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ABSTRACT 

We present chemical abundance measurements for 47 damped Ly-a systems (DLAs), 30 at z > 4, observed 
with the Echellette Spectrograph and Imager and the High Resolution Echelle Spectrometer on the Keck tele- 
scopes. H I column densities of the DLAs are measured with Voigt profile fits to the Ly-a profiles, and we 
find an increased number of false DLA identifications with SDSS at z > 4 due to the increased density of 
the Lyman-a forest. Ionic column densities are determined using the apparent optical depth method, and we 
combine our new metallicity measurements with 195 from previous surveys to determine the evolution of the 
cosmic metallicity of neutral gas. We find the metallicity of DLAs decreases with increasing redshift, improv- 
ing the significance of the trend and extending it to higher redshifts, with a linear fit of -0.22 ± 0.03 dex per 
unit redshift from z = 0.09-5.06. The metallicity "floor" of w 1/600 solar continues out to z ~ 5, despite our 
sensitivity for finding DLAs with much lower metallicities. However, this floor is not statistically different 
from a steep tail to the distribution. We also find that the intrinsic scatter of metallicity among DLAs of ~ 0.5 
dex continues out to z ~ 5. In addition, the metallicity distribution and the a/Fe ratios of z > 2 DLAs are 
consistent with being drawn from the same parent population with those of halo stars. It is therefore possible 
that the halo stars in the Milky Way formed out of gas that commonly exhibits DLA absorption at z > 2. 

Subject headings: galaxies: abundances — galaxies: ISM — galaxies: evolution — galaxies: general — 
Galaxy: halo — quasars: absorption lines 



1. INTRODUCTION 

The recent detection of star-forming g alaxies out to z ~ 8 
dBouwens et al.l2010l : lBunker et al.l2()10) provides new infor- 
mation about the crucial formative stages of galaxy evolution. 
One of the important goals of these studies is to identify the 
high-redshift progenitors of current galaxies and their stellar 
populations. A valuable diagnostic tool that links the prop- 
erties of high-redshi ft galaxies to their descendants is their 
chemical properties (Pettini 2006). While chemical proper- 
ties are exceedingly difficult to extract from the faint starlight 
emitted by star-forming galaxies at high redshift, they are 
more easily determined from protogalactic gas det ected in ab- 
sorption against bright background quasars (e.g. Pettin l2006l ; 
iProchaska et al.|[2"003al) . If the metallicity of such gas were to 
resemble the metallicity of, say, globular clusters, then a rea- 
sonable first guess is that stars in globular clusters form out of 
this gas. The evolution of the metallicity and other properties 
of this gas could then provide information about the sequence 
of events leading to the formation of current galaxies. 

The high-redshift Universe is occupied by populations of 
objects spanning a wide range of metallicities. While the 
mean metal abundances of dense , compact objects such a s 
quasars is given by [M/Hfl ~ (Ha mann & FerlandH l999). 
diffuse large-scale confi gurations such as t he Ly-a forest 
clouds have [M/H]=-3 dAguirre et all 120041) . In fact, ob- 
servers have no w even identified gas without any trace of 
heavy elements (Fumagalli et al. 2 01 lal) . With a mean metal 
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abundance of [M/H]=-1.4 dPettinil 120061). the dam ped Ly-a 
absorption systems (DLAs; see [Wolfe etallllOOll for a re- 
view) are galactic-scale objects that lie between these two ex- 
tremes. The purpose of this paper is to present the results of 
a new survey for the metal abundances of DLAs out to z ~ 5. 
We focus on DLAs for several reasons. First, since they are 
detected in absorption, DLAs are unbiassed with respect to 
luminosity and presumably mass. Second, in contrast to all 
other classes of quasar absorption systems, the gas in DLAs 
is mainly neutral. In fact, DLAs do minate the neutral-gas 
content of the Unive rse out to z ~ 5, (IProchaska et al.ll2.005t 
lO'Meara et al.l 120071) . and exhibit properties indicating that 
they are neutral-gas reservoir s for star formation at h igh red- 
shift (|Nag]unine2LiD|2^ More- 
over, the large optical depth at the Lyman limit of the neutral 
gas in DLAs (typically tll ~ 10 4 ) eliminates the need for un- 
certain ionization corrections to deduce the metal abundances. 
Accordingly, the most acc urately determine d abundances at 
high redshift are for DLAs dWolfe et al.12 005). Third, because 
DLA ga s likely originates in the outer regions of high-z galax- 
ies (e.g. Rafelski et al. 2011), chemical and dynamical infor- 
mation would be preserved due t o the long time-scales for dis- 
sipati on and dynamical mixing dFreeman & Bland-Hawthornl 
12001 . 

Our primary goal is to measure the cos mic metallicity, 
(Z) = logio(r2 m etai s /^gas) - log lo (M/H) dLanzetta et all 
1 1995b . at redshifts z > 4. Previously we determined (Z) by 
measuring the column-density wei ghted mean metal abun- 
dances (see Prochas ka et al.l l2003b) in redshift bins contain- 
ing equal numbers of DLAs. The data revealed a statistically 
significant decline in (Z) with increasing redshift, given by 
(Z) = (-0.26±0.07)z-(0.59±0.18), w hich was confirmed at 
z < 4 by i ndependent measurements (Kul karni et al.l 120051 
2007, 2010). While the evolution of (Z) appears linear in 
redshift, it is non-linear in time, with the slope, d(Z)/dt, be- 
coming significantly steeper at earlier cosmic times. Since the 
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iProchaska et al.l(l2003bl) sample contained only 10 DLAs with 
z > 4, (Z), and its derivatives were poorly determined at such 
high redshifts. Therefore, we aim to determine whethe r this 
trend continues back to z ~ 5. IProchaska et al.1 (l2003bl) also 
reported (1) a "metallicity floor" at [M/H]=-3.0, because no 
metal abundances were found below this value despite their 
sensitivity to [M/H] values well below -4.0, and (2) a large 
dispersion of [M/H] at all redshifts. We wish to investigate 
whether these phenomena persist at higher redshift s, as at 
some p oint in time the gas should be pristine (Fumagalli et al. 
1201 lah . Although increasing the redshift range of our survey 
to z ~ 5 corresponds to an increase in survey time interval of 
only Ri 0.5 Gyr, we shall show that this is a critical time inter- 
val in which, as stated above, d(Z)/dt becomes increasingly 
steep with decreasing t. 

In addition, we wish to assess if there is a systematic bias 
in the Sloan Digital Sky Survey (SPSS) N m statistic s at z > 4 
(IProchaska & Wolfd 2009: No terdaeme et al.ll2009h . such as 
the column density distribution function, total covering frac- 
tion, and cosmological Nm mass density. These statistics are 
determined solely from the discovery spectra of the SDSS 
database, with low resolution of FWHM^ 2A and sometimes 
at low signal-to-noise (S/N) (with a cut at S/N> 4). The anal- 
ysis of the damped Ly-a profile is significantly affected by 
line-blending with the Ly-a forest 'clouds', and the increase 
of the density of the Ly-a forest at high redshift may lead 
to an overestimate of the Nni measurements and an increased 
false positive identification rate fo r DLAs. While such biases 
are likely to be small at 2 < z < 4 (Prochaska et al. 2005), our 
new measurements resolve out the Ly-a forest and enable the 
constraint of such a bias at z > 4. 

Another goal is to compare the metallicity distribution of 
DLAs with those of known stellar population s. Previous 
studies dPettini et al.l 1 1997h IPettinil 120041 12006ft have shown 
that the distribution of metallicities at z > 2 disagree with 
known stellar populations in the Galaxy. However, if high 
redshift galaxies are surrounded by a thick disk of neu- 
tral gas (IProchaska & Wolfelll997bl:IWolfe etal.ll200l 120081: 
iRafelski et al.ll2009l.l201 lblFumagalli et al.ll201 lbl) . then stars 
that formed in situ in such gas at high redshift should be 
present in z = galaxies, such as the Milk Way. These stars 
would presumably have a similar metallicity distribution as 
the DLA gas. We wish to investigate whether or not this dis- 
agreement persists in our sample that is larger in size and ex- 
tends to higher redshifts. 

The outline of this paper is as follows: The observations 
and data reduction are presented in §|2] and the analysis of the 
resultant spectra is detailed in §|3] In §13.11 Voigt profiles are 
fitted to the damped Ly-a profiles to determine the H I col- 
umn densities, A^hi- We discuss the possible elemental abun- 
dances that can be used to determine the metallicity of DLAs 
in §13.21 Column densities of the metal l ines are determined 
from t he apparent optical dep th method (ISavage & Sem bach 
I199U IProchaska eta l. 2001b), and combined with the N m 
values to determine element abundances for the DLAs in §13.31 
This sample is augmented by unbia sed a bundances obtained 
from the literature as described in § 13.41 We justify the use 
of vario us elements for determining the metallicities of DLAs 
in §13.51 The determ inati on of (Z) as function of redshift and 
time is presented in §13.61 and a comparison of the DLA metal- 
licity distribution with that of halo stars is shown in §13.71 In 
§H] we discuss the metallicity distribution and evolution, and 
its implications for halo star formation. Specifically, we com- 



pare recent si mula tions to our data and discuss the metallicity 
evolution in §14.11 and consider the implications of the metal- 
licities and a-enhancement of DLAs for halo star formation 
in §14.21 In §J5]we summarize the results. 

Throughout this paper we adopt a cosmology with 
(fi M) ft A ,/0=(0.3,0.7.0.7). 

2. OBSERVATIONS AND DATA REDUCTION 

The primary scientific goals of our survey are to: (i) per- 
form an unbiased evaluation of the metallicity in neutral 
atomic-dominated hydrogen gas at z > 4; (ii) confirm and re- 
fine previous estimates on the incidence and H I mass density 
in z > 4 DLAs; and (iii) identify new targets for follow-up 
ALMA observations in the C II* 158 micron transition. These 
goals were well suited to the spectroscopic cha racteristics of 
the E chellette Spectrograph and Imager (ESI; Sheinis et al. 
2002), with the possibility of followup observatio ns with the 
High Resolution Echelle Spectrometer (HIRES; IVogt et al] 
1994) on the Keck I and II 10m telescopes. We therefore re- 
quired a set of z > 4 DLAs that could be observed from Keck 
observatory. 

The first sample of targets were drawn from the SDSS-DR3 
an d SDSS-DR5 surveys fo r da mped Ly-q systems perform ed 
bv lProchaska et al.1 (|2005) and IProchaska & Wolfd (120091) re- 
spectively. We restricted the list to sources that were suffi- 
ciently bright such that we could obtain S/N > 1 5 spectra in a 
non-prohibitive observing time. After the first few nights, we 
needed to supplement the sample with additional z > 4 DLAs, 
and we searched for new z > 4 quasars from SDSS-DR7. 
These additional targets will not be included in the H I analy- 
sis as they were selected in a more heterogeneous manner than 
the original set. They are unbiased, however, with respect to 
the gas metallicity because they were selected solely on the 
presence of strong H I Ly-a absorption. In addition to this 
primary high redshift quasar sample obtained in 2009-2011, 
we also include new observations of lower redshift quasars 
(2.5 < z < 3.5) obtained in 2008 to augment the number of 
unbiased metallicity measurements (see below). We note that 
the effective upper limit to the redshifts of our DLA sample 
is about 5.3; DLAs with higher redshifts are essentially unde- 
tectable because line blanketing by the Ly-a forest severely 
attenuates the quasar continuum at wavelengths shortward of 
Ly-a emis sion; i.e. at wavelengths where such DLA s would 
be located dFan et al.ll2006l:ISongaila & Cowidl2010l) 

The new observations comprise a total of 6 nights (effec- 
tively 5 nights when considering the losses to poor weather) 
on ESI and 7 nights on HIRES. In Table Q] we present a jour- 
nal of observations describing the 50 quasars in our sample. 
Column (1) gives the quasar coordinate name obtained from 
the SDSS survey. Columns (2) and (3) give the RA and Dec. 
of the quasar in J2000. Columns (4) and (5) give the AB r 
and i magnitudes of the quasar: in most cases the i magnitude 
gives a realistic measure of the unattenuated quasar bright- 
ness since at the redshifts of most of the DLAs, the i filter is 
placed redward of Ly-a emission. Column (6) gives the emis- 
sion redshift of the quasars, and column (7) lists candidate 
DLA absorption redshifts obtained from SDSS. Column (8) 
provides the UT date of the observations, column (9) states 
the instrument used, and column (10) lists the slit size. Col- 
umn (11) gives the exposure time, and column (12) the S/N 
pixel" 1 of the data. 

In general, we obtained 2-3 exposures per quasar with ESI, 
with sufficient exposure time to reach S/N values larger than 
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Table 1 

Journal of Observations 





K.A. 


Dec. 


r 


i 


Zem 


Zabs 


Date 


Instr. 


Mil 


U'.xp 


o/JN 




(J 2000) 


(J 2000) 


(AB) 


(AB) 






(UT) 




(arc sec) 


(s) 


pixel -1 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


J0040-0915 


00 40 54.7 


-09 15 27 


20.48 


19.18 


4.98 


4.74 


2010 Jan 20 


ESI 


0.50 


2280 


16 
















2011 Jan 16 


HIRESr 


0.86 


7200 


6 
















201 1 Jan 24 


HIRESr 


0.86 


3600 




J0210-0018 


02 10 43.2 


-00 18 18 


20.48 


19.20 


4.72 


4.57 b 


2010 Jan 20 


ESI 


0.50 


2880 


9 


J033 1-0741 


03 31 19.7 


-07 41 43 


20.56 


19.12 


4.71 


4.19 b 


2010 Jan 20 


ESI 


0.75 


4320 


18 


J0747+4434 


07 47 49.7 


+44 34 17 


19.62 


17.35 


4.43 


4.02 


2009 Mar 22 


ESI 


0.75 


3600 


15 
















2011 Jan 16 


HIRESr 


0.86 


3601 


9 
















201 1 Jan 24 


HIRESr 


0.86 


10800 




J0759+1800 


07 59 07.6 


+18 00 55 


21.06 


19.16 


4.86 


4.66 


2010 Jan 20 


ESI 


0.75 


4320 


18 


J08 13+3508 


08 13 33.3 


+35 08 11 


20.81 


19.13 


4.92 


4.30 b 


2010 Apr 21 


ESI 


0.75 


1400 


14 



Note. — Units of right ascension are in hours, minutes, and seconds, and units of declination are in degrees, arcminutes, and arcseconds. Photometry 
from SDSS l Abazajian et al. 2009). This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for 
guidance regarding its form and content. 



a Median signal-to-noise per pixel at A S3 8300A, except in the HIRES observations of J1353+5328 and J1541+3153, where it is at A S3 7300 A. 
b Marked as a DLA in SDSS, but is not a DLA based on higher resolution ESI data 



15 pixel" 1 (11 km s" 1 ). Observations were reduced on the fly 
to determine if the DLA candidate did not meet our thresh- 
old criteria for being a DLA, i.e. N(HI)> 2 x 10 20 cm" 2 . If it 
did not, then we abandoned obtaining more exposures, result- 
ing in reduced S/N values for those targets. These quasars are 
nonetheless listed in TableQ] and the corresponding DLAs are 
marked with the footnote b. TableQ]also lists the median S/N 
per pixel at A « 83OOA0, corresponding to the typical wave- 
length of the measured metal absorption lines. The average 
S/N of the ESI spectra with confirmed DLAs is 25, with a 
minimum S/N of 14. The mean S/N of the spectra with DLA 
candidates not meeting our criteria is 15, with a minimum S/N 
of 9. 

All ESI observations prior to 2010, and all HIRES obser- 
vations were observed at the same slit position to facilitate 
the removal of cosmic rays between exposures. In order to 
decrease the effects of fringing in ESI observations (~ 15% 
at 9000A; ISheinis et ail 120021) at the longest wavelengths 
(A > 8,000A), we used a three point dither pattern along the 
slit for all ESI observations starting in 2010. Throughout the 
ESI observations, we primarily used the 0."75 slit, except in 
two cases when we used the 0."5 slit, which correspond to a 
resolution of FWHM of ~ 44 and ~ 33 km s" 1 respectively 
(R « 7,000-10,000). For the HIRES observations, we used 
the CI decker, with a 0."86 slit corresponding to a resolution 
of FWHM of 6.2 km s" 1 or 1.5 km s" 1 pix" 1 (R « 50,000). 
ESI has a fixed format echellette that covers the spectral re- 
gion A = 4,000 to 11,000A, while the HIRES echellette cov- 
erage varies with each setup. The HIRES observations are all 
followup observations of the ESI observations in order to de- 
blend lines, check for saturation, and measure weak lines, and 
we therefore targeted the appropriate metal lines for metallic- 
ity determination with the HIRES setups. 

Figure Q] shows the 25 quas ar sp ectra containing 32 con- 
firmed DLAs with z > 4 (see 83. II regarding confirmation of 
DLAs), in order of increasing value of the largest redshift 
DLA for each spectrum. We note that while these spectra also 
contain other lower redshift DLAs, we only show those with 
z > 4. In addition to the 25 quasars containing z > 4 DLAs, 

5 The HIRES observations of J1353+5328 and J1541+3153 have their S/N 
determined at A ss 7300A due to their lower redshifts and corresponding 
HIRES coverage. 



we also analyze ESI observations and HIRES followup from 
2008 which were obtained in a similar unbiased fashion. 

All the data were reduced with the XIDI0 package devel- 
oped in IDL by J. X. Prochaska. Spec ifically, the ESI observa- 
tions were reduced with ESIRedtnQ (Prochaska e t al.ll2003ah 
and the HIRES observations were reduced with HIReduxQ, 
both of which are publicly available. All the data were ex- 
tracted with optimal extraction, and the data were continuum 
fit using the x_contimmm routine within XIDL. 

3. ANALYSIS AND RESULTS 

In this section we present the analysis of the quasar spec- 
tra obtained in this survey. We first confirm the DLA candi- 
dates using Voigt profile fitting, determine their TVhi values, 
and investigate the false positive rate of SDSS-selected DLAs 
at high redshift. We discuss the possible elemental abun- 
dances that can be used to determine the metallicity of DLAs, 
and then describe the derivation of elemental abundances. We 
then compile the latest literature metallicities, justify the use 
of multiple elements for metallicity determination, and deter- 
mine (Z) as function of redshift. Finally, we compare the 
metallicity distribution of DLAs to that of known stellar pop- 
ulations in the Milk Way, focusing on the halo stars. 

3.1. HI Column Densities 

We determine the Nui values of the damp ed Ly-a systems 
in our s ample using the same methodology as Prochaska et al. 
(2003a). Specifically, we fit the Ly-a line with a Voigt profile 
in the fluxed, but not continuum fit, spectra using the routine 
x_f itdla, part of the XIDL package. We tie the centroid of 
the Voigt profile to the redshift of the low-ion metal-line tran- 
sitional which yields a velocity centroid to within ~ 30 km 
s -1 of the true line center (Prochaska et al. 2003b). We simul- 
taneously fit the Voigt profile and the continuum to obtain a fit 
that matches both the core and wings of the profile. The un- 
certainties associated with the redshift errors are significantly 
less than the uncertainty in the continuum fit, and therefore are 
not separately included in the fitting procedure. While this fit- 
ting procedure is not completely quantitative, it is the standard 

6 See |http : / /www . ucolick . org/~xavier / IDL/| 

7 See ht tp : / / www2 . keck . hawaii . edu/ inst /esi/ESIRedux/| 

8 See |h^^^^^^^^^^K^^^l^^^^^^^ca^r^e^^^^e'dhj^^| 

9 All of the bona-fide DLAs exhibit at least one positive detection. 
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Figure 1. ESI spectra of 25 quasars with 30 DLAs at z > 4. Fluxes (black curves) and 1-tr error arrays (green curves) are normalized to unit continuum. The 
orange lines are zero flux, z > 4 Damped hya lines, denoted by red arrows, are listed in order of increasing value of the largest redshift DLA for each spectrum 
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methodology and is justified by the fact that the uncertain- 
ties are dominated by systematic errors: continuum fitting and 
line-blending. In our experience, a st andard \ 2 analysis yield s 
unrealistically small error estimates dProchaska et al.l l2003a). 
We place very conservative error estimates on the Nm val- 
ues by manually selecting values that ensure that the error re- 
gion encompasses all possible solutions for A^hi- These values 
systematically increase with increasing redshift due to the in- 
creased Ly-a forest blending, and the reader should consider 
them to be roughly at the 95% confidence level. 

The sample of candidate DLAs all have A^hi> 2 x 10 2l) cm~ 2 
based on the SDSS spectra with spectral resolution FWHM 
— 2A. At low redshift, this is sufficiently high resolution to 
obtain relatively accurate Nm, and therefore yields only a 
small number of false candidates. However, as we move to 
higher redshift, the Ly-a forest becomes significantly more 
dense, and the damped Ly-a profiles in the SDSS spectra are 
significantly affected by line-blending with random Ly-a for- 
est clouds. Since the ESI spectra resolve out much of the Ly-a 
forest, we can therefore obtain more robust values of Nhi- In 
order to confirm a candidate DLA from SDSS to be a true 
DLA, we require Nhi> 2 x 10 2() cm" 2 in the ESI spectrum. 

We fitted Nm Voigt profiles to 68 candidate DLAs from 
SDSS, and confirm 51 as bona-fide DLAs. Figures |2] and [3] 
show the Voigt profile fits of 5 1 newly confirmed DLAs or- 
dered by increasing RA. The velocity centroid of the DLA 
is determined from the low-ion metal transitions, and some 
panels have more than one DLA candidate marked, and in 
those cases the central DLA at zero velocity is the DLA for 
the information given in that panel. The velocity interval of 
the individual panels varies with column density such that the 
wings of the damped profile are presented. 

The higher density of the Ly-a forest at higher redshifts 
results in more line blending than at lower redshift, and there- 
fore sometimes one or both of the wings are not clearly de- 
fined, and in these cases we rely on the core of the line. In a 
few other cases, a strong Ly-a absorber is blended with one 
side of the DLA's core, and the resultant profile is constrained 
by the core on the other side, and the wings on one or both 
sides. For example, the Voigt profile fit of J 1200+46 18 at 
Zabs = 4.477 turns up more sharply than the data on the blue 
side of line center due to confusion with one or more Ly-a 
absorbers at v ~ -900 km s" 1 . The same effect can be seen in 
J1201+21 17 at z abs = 4.158, and a few other DLA fits. These 
challenges are not a problem, as the redshift is determined 
from the low-ion metal lines, and we therefore have sufficient 
information for a robust determination of A^hi- 

The 17 false candidates consist of multiple lower column 
density lines that are blended together in the lower resolu- 
tion SDSS data. Of the 68 candidate DLAs in Table [TJ 46 
of them have z > 4. Of those 46 candidate z > 4 DLAs, 14 
are false positives, leaving a total of 32 confirmed DLAs with 
z > 4. As a result, we find that —30% of z > 4 DLAs are 
misidentified in SDSS due to line-blending. In contrast, for 
our sample of 22 z < 4 DLAs, we find 3 false candidates, 
none of which are one of the 14 z ~ 3 candidate DLAs ob- 
tained in our 2008 run. To investigate this further, we limit 
ourselves to the more homogeneous part of the sample that 
were ident ified with the a utomated routine from DR5 SDSS 
data dProchaska & W olfe 2009), excluding the manually se- 
lected DR7 SDSS DLAs (see C3>. Also, we exclude 2 DLAs 
that were lower than N m > 2 x 10 2() cm" 2 in the SDSS data, but 
were measured in the ESI spectra to meet the DLA criterion. 



This removes 5 false candidate DLAs (4 of which are at z > 4) 
and 9 confirmed DLAs from the sample, leaving a total of 56 
candidate DLAs. Out of these 38 candidate z > 4 DLAs, 10 
are false positives, leaving a total of 28 confirmed DLAs with 
z > 4, and a —26% misidentification rate. For the 18 z < 4 
DLAs, we find 2 false candidates, and a —11% misidentifi- 
cation rate. We note that one of these two DLAs is found in 
a quasar with z em = 4.9, and therefore is severely affected by 
the z > 4 Ly— /3 forest, which is uncommon, and therefore the 
misidentification rate at z < 4 is likely closer to —6%. Despite 
the small number statistics, there clearly is a larger percentage 
of false positives at z > 4 than at z < 4. 

We compare the Nm values determined from the DR5 SDSS 
data and from our ESI data in Figure |4] The top panel is 
the ESI Nm versus the SDSS Nm, while the bottom panel 
are number counts of misidentified DLAs. The dotted line 
marks the minimum criteria to be a DLA, A^hi> 2 x 10 20 cm" 2 . 
There is general good agreement between the ESI Nm values 
and the SDSS Nm values within their uncertainties. In addi- 
tion, the histogram shows that the primary misidentifications 
occur for DLAs near the Nm= 2 x 10 20 cm" 2 limit. Some of 
these misidentified systems are super Lyman limit systems, 
although the majority are due to a combination of lines in 
the Lyman-a fo rest. Studies using DLA id entifications from 
SDSS at z > 4 dProchaska & Wolfdl2009h therefore need to 
apply a correction for the increased false positive rate before 
determining the column density distribution function, the to- 
tal covering fraction, and the integrated mass density of H I 
gas. The histogram in Figure |4] suggests that this correction 
mainly needs to be applied for the lowest Nm DLAs. 

For such future corrections, we break down the —26% 
misidentification rate into random and systematic misidenti- 
fications. Of the 10 false positives from DR5, 3 are random 
misidentifications with high Nm absorbers at the correct red- 
shift that do not quite meet the criteria of A^hi> 2 x 10 20 cm" 2 , 
but are within 1 - 2a of their previous measurements. The 
other 7 false positives are systematic misidentifications due to 
a combination of many lower Nm absorbers masquerading as 
a single higher A^hi absorber. Figure |5]show Voigt profile fits 
to both the original SDSS data and the higher resolution ESI 
data for the 7 systematic misidentifications from DR5. 

We expect that a number of absorbers with Nm< 2 x 
10 20 cm" 2 will randomly be measured above this value, and 
that a number of absorbers with Nm> 2 x 10 20 cm~ 2 to be 
measured below this value. We check if the number of ran- 
dom misidentifications measured is as expected from the un- 
certainties and the column densi ty distribution function using 
the E ddington bias formalism ( Edding ton! 1 19 131: iTeer ikorpi 
120041) . We limit ourselves to DLAs with Af H i< 20.5, as above 
this Nm the problem is small (see Figure |4]i and a random er- 
ror is unlikely to move an absorber with Nm< 2 x 10 2() cm~ 2 to 
a higher Nm- The observed number of DR5 DLAs with z > 4 
and Nm< 20.5 is 9, to which we add the 3 random misiden- 
tifications for a total of 12 absorbers. We use an uncertainty 
of 0.2 dex, and model the column density distribution func- 
tion. The expected number of DLAs from the Eddington bias 
formalism is 8.4, resulting in an observed difference of 3.6 
DLAs. This is consistent with the 3 random misidentifica- 
tions measured, especially given the small number statistics 
considered here. 

We note that in our survey for DLAs, we reject proxi- 
mate DLA systems that are found to be within 3000 km 
s" 1 of the quasar emission redshift, as these DLAs may be 
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Figure 2. Voigt profile fits to 51 confirmed DLAs ordered by RA, with the x-range based on the iVm. Blue curves are best-fit profiles and gray shade includes 
95% confidence limits surrounding best fits. The red dotted line marks the velocity centroid of the DLA determined from the low-ion metal transitions, with 
possible multiple components as determined by the metal lines. Individual panels include quasar name, DLA redshift and H I column density. The green line 
represents the uncertainty. 
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Figure 4. Comparison of the Nya values determined from the DR5 SDSS 
data and from our ESI data in Figure [2] The top panel is the ESI N\\\ vs. 
the SDSS Nm, while the bottom panel are number counts of misidentified 
DLAs, based on the higher resolution and higher S/N ESI data. The dotted 
line marks the minimum criterion to be a DLA, N\\\> 1 X 10 20 citT 2 . 

phy sically associated or otherwise affected by the quasar 
(see IProchaska et al.ll2008t lEllison et ail 12010. EoTTl) . This 
removes our highest redshift confirmed DLA (z = 5.18, 
J1626+275), and it is not included in the analysis of the metal- 
licity evolution, although we include it for completeness in 
our tables. 

3.2. Determining DLA Metallicities 

A principal goal of this paper is to trace the build-up of 
heavy elements in neutral gas across cosmic time. In the Sun, 
oxygen (O) dominates by numbe r and contribut es w 50% of 
the heavy element mass budget (Lodders 2010). This frac- 
tion is likely even higher in the early universe when Type la 
supernovae h ave had less time to enrich the iron (Fe) peak el- 
ements (see 33.71 ). Indeed, O is commonly observed to be en- 
hanc ed in metal-poor en vironments such as our Galactic halo 
(e.g. lBesselletal][T99l . Ideally, one would assess the heavy 
element mass density through measurements of O. Unfortu- 
nately in DLA analysis, the O I transitions are generally ei- 
ther too weak or too strong to permit a precise estimate of its 
abundance relative to hydrogen. After oxygen, carbon (C) and 
neon (Ne) contribute ps 30% of the heavy element mass in the 
Sun. The same challenges as for oxygen apply to C, and there 
are no Ne transitions available for DLAs. It is evident that we 
are forced to use a proxy for the heavy element abundance. 

The elements nitrogen (N), magnesium (Mg), silicon (Si), 
sulfur (S), and Fe each contribute ks 5% of the heavy element 
mass density in the Sun and each of these can, in principle, 
be measured in DLAs. To best trace the dominant elements 
(O, C, Ne), one would prefer an element that is produced 
and expelled in tandem, i.e. predominantly in the supernovae 
explosions of massive stars. Both Si and S are good can- 
didates since they lie along the a-element sequence of nu- 
cleosynthesis, are produced in massive stars, and are easily 
measurable in high redshift DLAs. S is a volatile element 
with negligible depletion and with multiple lines of varying 
oscillator strengths. Alas, these lines are often in the Ly-a 
forest, making them measurable in only a fraction of DLAs. 
Si has many lines of varying oscillator strengths at multiple 
wavelengths, which results in at least one line being unsatu- 
rated and outside the forest for the majority of DLAs. While 
Si has nearly twice the mass density of S, it is a mildly re- 
fractory element and its gas-phase measurement may under- 



estimate the total m etal abundance (Prochaska & Wolf9 l2002T: 
IVladilo et al.ll201 lb . Considering all of these aspects together, 
we conclude that S is the best option and use its abundance as 
a proxy for the heavy element mass density in DLAs wher- 
ever possiblfU Under the expectation that metal-poor gas is 
relatively undepleted, we adopt Si as the second choice, espe- 
cially in low metallicity system s. W e further investigate the 
possibility of depletion of Si in 33.51 where we determine that 
it is not significantly depleted in DLAs based on current mea- 
surements. 

Within the DLA literature, zinc (Zn) has played a spe- 
cial role in estimating the metallicity (e.g. | Mever et al1ll989t 
iPettini et alJ fl990l fl99i \l99% lAkerman et al l 120051) . This 
reflects both observational convenience and the fact that Zn 
is non-refractory. Specifically, Zn has two strong transitions 
at 2026 and 2062 A, that almost always lie outside the Ly- 
a forest, and are rarely saturated due to their low oscillator 
strengths and low Zn abundances. Unfortunately, Zn is only 
a trace element which contributes w 10~ 4 of the mass den- 
sity for the heavy elements. Furthermore, the weakness of 
the Zn transitions make it impossible to measure its abun- 
dance in low metallicity DLAs, and the large rest-frame wave- 
lengths of the transitions make them difficult to measure at 
high redshift. In addition, the nature of Zn is not fully un- 
derstood; Zn is not an a-element, nor an Fe peak element. 
While Zn is often observed to track Fe, theoretical work on 
the nucleosyn thesis of Zn often predic t s it may originate in 
massive stars dWooslev & Wea ver 1995; Hoffman et al.ll996t 
Umeda & Nomoto 2002; iFenner et all 12004). On the more 
positive side, Zn tracks the relative abundances of S and Si in 
DLA s, an d empirically is a fair substitute for these elements 
(see 33.5b . Therefore, when O, S, and Si are either not cov- 
ered by our spectra or are saturated, and the Zn transitions are 
available, we use the Zn abundance as a proxy for the heavy 
element mass density. 

In the small number of cases where none of the above el- 
ements are available, we are left with using elements such 
as Fe, nickel (Ni), or aluminum (Al). As Fe and Ni are re- 
fractory elements, they will be depleted at high metallicities. 
At low metallicities, we can measure the intrinsic abundance 
of Fe and Ni, which tracks the abundances of S and Si well, 
with an offset due to the a-enhancement of DLAs. We deter- 
mine the a-enhancement in 33.3l and 33.5l from the abundance 
of Fe combined with an a-element abundance. We apply an 
[a/Fe] offset of 0.3 dex, which is slightl y different than that 
used i n some previous studie s (0.4; e.g. IProchaska & Wolfd 
l2002t [Prochaska et al. 2003b), but i s more representativ e of 
[a/Fel of DLAs dWolfe etall 120051 and see below in 3141 
and {230 In addition, we also determine [Fe/H] from iron 
peak elements for all DLAs, as this is useful for studying the 
chemistry of DLAs, and for an estimate of the dust-to-gas ra- 
tio, as Fe is a dep leted element for DLAs with [M/H] > -1 
(Wolfe etal. 200J). In cases where Fe is not available, we de- 
termine [Fe/H] useing Ni with a 0. 1 dex correction, or Al. We 
note that whenever we use an Fe peak element for the metal- 
licity, we add the uncertainty in [a/Fe] of 0.16 to the [Fe/H] 
uncertainty in quadrature. 



Aside from the few cases where oxygen is measured precisely. 
1 1 We note that the y-axis of Figure 8 in Wolfe et al. 1 2005) is mislabeled. 
The tick labeled 0.5 should be labeled 0.7 
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3.3. Measured Element Abundances 

In this subsection we describe the derivation of ionic col- 
umn densities for all of the new confir med DLAs using the 
appar ent optical depth method (AODM; ISavage & Sembach] 
1 1 99 11) . This method can uncover hidden saturated lines by 
comparing the apparent column density per unit velocity, 
N a (v), for multiple transitions of a single ion. We then cal- 
culate N a (v) for each pixel from the optical depth equation, 
and sum over the velocity profile of each transition to ob- 
tain the total column density, N. If multiple lines for a sin- 
gle ion are available, then we take the weighted mean of the 
column densities from unsaturated lines. The AODM tech- 
nique yields column densities that agree well with line fitting, 
and it is easier to apply to a large data sets such as this one 
dProchaska et a l. 2001b). We use the wavelength s and oscilla- 
tor strengths shown in Table 2 of lMortonl (12003|) and the me- 
teoritic solar abundances shown in Table 1 of lAsplund et al] 
(2009). 

The medium resolution of ESI could result in possible satu- 
ration of unresolved weak metal lines if they were particularly 



narrow (Penprase et al. 2010), and such lines should normally 
be treated as lower limits. However, in this study we obtained 
HIRES followup in order to 1) check for saturation of nar- 
row weak metal lines from the ESI measurements by resolv- 
ing the lines, and 2) obtain metallicities that would otherwise 
not be measurable with ESI due to blending of lines or weak 
transitions. By targeting the most likely problematic candi- 
dates with HIRES, we reduce the saturation issues. Specifi- 
cally, we followed up DLAs with narrow lines whose metal- 
licity was based on single metal lines, as these are the most 
likely culprits. In the cases where a line looked potentially 
saturated without available HIRES data, we used a different 
line instead. In a few cases, this resulted in an Fe-element 
based metallicity rather than an a-element based metallicity. 
We therefore use the remaining ESI sample without apply- 
ing saturation corrections. In addition, if we did miss a small 
number of saturated DLAs, they would not affect the main 
conclusions, as we have sufficient measurements per redshift 
bin such that (Z) is not sensitive to any individual measure- 
ment (see §3.4). Lastly, as we probe higher redshifts, the 
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Figure 6. DLA metal abundance, [M/H], versus redshift for the 46 DLAs in 
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Figure 7. DLA metal abundance, [M/H], versus redshift for 195 DLAs from 
the literature. The symbols and colors indicate the origin of the [M/H] values: 
Si II measurements (green filled circles), S II measurements (blue triangles), 
Zn II measurements (brown boxes), Fe II measurements (red crosses), and 
others such as O I (gold stars) and the combination of limits (black diamonds). 



common element transitions used for metallicity determina- 
tion in DLAs are less likely to b e saturated due to their l ower 
metallicities at higher redshifts (Proch aska et al.ll2003bl and 
this study). 

In Table|2]we present new metallicity measurements for 47 
DLAs, with 30 of them at z > 4. Four DLAs do not have mea- 
surable metallicities due to a lack of unsaturated relevant lines 
red ward of the Lyman-a forest, and are not in the table. We 
do note, however, that each of these exhibits obvious metal ab- 
sorption and therefore has a metallicity [M/H] > -2. Another 
is a proximate DLA, and is included in the table. Column 
(1) gives the quasar coordinate name obtained from the SDSS 
survey, column (2) gives the DLA absorption redshifts deter- 
mined from the metal transition lines, and column (3) gives 
the logarithm of the measured Nhi values. Columns (4) and 

(5) give the alpha element flag and metallicity, and column 

(6) and (7) give the Fe flag and metallicity. Column (8) and 
(9) give the flag and final metallicity used in this paper. The 
flags in columns (4) and (6) state which transitions are used 
and if they are detections or limits, as described in the Table 
notes. We also provide the measured elemental column den- 
sities of these DLAs in the online material in the appendix. 
In addition, we show the velocity profiles for those metal-line 
transitions in the online material. 

The uncertainty in [o/H], [Fe/H], and [M/H] in Table[2]in- 
clude the uncertainty of A^hi added in quadrature, and there- 
fore are dominated by the uncertainty in A^hi- We ignore the 
error in the continuum fit in our error calculations for metal - 
line transitions, although it may be important for very weak 
transitions. We set a minimum error of 0.1 dex for all metal- 
licity measurements, and note that our main results are not 
sensitive to the uncer taint ies in individual measurements (see 
bootstrap method in ^3. 6b . 

In Figure [6] we plot the 46 metallicities obtained in this 
paper as a function of redshift (not including the proximate 
DLA). The figure outlines which elements are used to ob- 
tain [M/H] and the uncertainties of those measurements. The 
symbols and colors indicate the origin of the [M/H] measure- 
ments, with blue triangles representing S II, green filled cir- 
cles representing Si II, and red plus crosses representing Fe II. 
We note that the metallicities derived from Fe II include an 
q-en hanc emen t correction, with [M/H] = [Fe/H]+0.3 dex (see 
Eland €DT >. 



3.4. Element Abundances from the Literature 

In this subsection we describe the compilation of metallic- 
ities of DLAs from the literature. Over a thousand DLAs 
have been identified by SDSS (e.g. [P rochaska et al.l 120051; 
iProchaska & Wolfel2009tlNoterdaeme et al.ll2009h .and many 
hundreds have been observed with a wide variety of tele- 
scopes and instruments with v arying resolution, S/N, and 
wavelength coverage. Similar to IProchaska et al.l (l2003bl) . we 
restrict our analysis to DLAs observed at high-resolution (R > 
5000) and high S/N (15 pixel" 1 ), and require them to meet the 
DLA criteria of Nm> 2 x 10 20 cm" 2 . In addition, we avoid 
samples with metallicity biases, such as metal poor sam- 
ples ( e.g. lPettini et al.1200 8: Penpra s"e et al.l20iotlCooke et al.l 
1 20111) or metal ric h samples (e.g. Herbert-Fort et al.l 120061 
Kapla neit alj I2010I) . We make an exception for low red- 
shift DLAs with z < 1.5, a s these are genera l ly found based 
on str ong metal lines (e.g. iKhare et al.ll2004t [Kulkarni et alj 
12005b iMeiring et ail 12006b IPeroux et al.l 120061) . Therefore, 
other than for DLAs in our lowest redshift bin, the literature 
sample is not (explicitly) metallicity biased. 

In Figure|7] we plot the compilation of 195 metallicities lit- 
erature abundance values as a function of redshift. The figure 
outlines which elements are used to obtain [M/H] and the un- 
certainties of those measurements. The symbols and colors 
indicate the origin of the [M/H] measurements, with blue tri- 
angles representing S II, green filled circles representing Si II, 
red plus crosses representing Fe II, brown boxes representing 
Zn II, and gold stars representing O I, and combination of lim- 
its (black diamonds). We also present these values in Table[3] 
where column (1) gives the quasar name, column (2) gives the 
DLA absorption redshifts, and column (3) gives the logarithm 
of the measured A^hi values. Columns (4) and (5) give the al- 
pha element flag and metallicity, and column (6) and (7) give 
the Fe flag and metallicity. Column (8) and (9) give the flag 
and final metallicity used in this paper, and column (10) cites 
the references. Similar to Table|3] the flags in columns (4) and 
(6) state which transitions are used and if they are detections 
or limits, as described in the Table notes. The uncertainties in 
the table include the uncertainty in A^hi added in quadrature. 

We combine the 195 literature abundance values with our 
new ESI and HIRES observations for a total of 241 abun- 
dances. We investigate the abundance distribution of DLAs, 
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Table 2 

New DLA metallicities 



QSO 




logNffl 




[a/K] 


/Fe" 


[Fe/H] 


/mtl c 


[M/H] 


\ 1 1 


\A> 




Id) 


(J) 




\ t ) 


(o) 


\ 7) 


J0040-0915 


4.7394 


20.3()±0.15 







1 


-1.70±0.16 


2 


-1.40±0.17 


J0747+4434 


4.0196 


20.95±0.15 


3 


<-1.90±0.15 


4 


-2.58±0.20 


2 


-2.28±0.20 


J0759+1800 


4.6577 


20.85±0.15 


4 


-1.74±0.16 


3 


<-1.14±0.15 


1 


-1.74±0.16 


J0817+1351 


4.2584 


21.30±0.15 


4 


-1.15±0.15 


1 


-1.30±0.16 


1 


-1.15±0.15 


J0825+3544 


3.2073 


20.30±0.10 


2 


>-1.83±0.10 


1 


-1.98±0.10 


2 


-1.68±0.16 


J0825+3544 


3.6567 


21.25±0.10 


1 


-1.83±0.13 


4 


-1.95±0.11 


1 


-1.83±0.13 


J0825+5127 


3.3180 


20.85±0.10 


1 


-1.67±0.14 


1 


-2.08±0.10 


1 


-1.67±0.14 


J083 1+4046 


4.3440 


20.75 ±0.15 


1 


-2.36±0.15 


1 


-2.41±0.17 


1 


-2.36±0.15 


J0834+2140 


3.7102 


20.85±0.10 


3 


<-1.69±0.10 


1 


-1.86±0.10 


2 


-1.56±0.16 


J0834+2140 


4.3900 


21.00±0.20 


4 


-1.30±0.20 


1 


-1.69±0.20 


1 


-1.30±0.20 



Note. — Note that none of the reported limits take into account the uncertainty in Nm- This table is available in its entirety in 
a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content. 



a = No measurement; 1= Si measurement; 2 = Si lower limit; 3 = Si upper limit; 4 = S measurement, . 

b = No measurement; 1= Fe measurement; 2 = Fe lower limit; 3 = Fe upper limit; 4 = [Ni/H]-0. 1 dex; 5 = [Cr/H]-0.2 dex; 6 = 
[Al/H]. 

c 1 = [a/H]; 2=[Fe/H]. In the latter case, we use [M/H] = [Fe/H]+0.3 dex. 

^ This is a proximate DLA, and is therefore not included in the analysis of the metallicity evolution. 



Table 3 

Literature DLA metallicities 



QSO 




log/V H i 




[a/H] 




[Fe/H] 


/mtl C 


[M/H] 


Refs d 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


Q2359-02 


2.0951 


20.70±0.10 


1 


-0.72±0.10 


1 


-1.61±0.10 


1 


-0.72±0.10 


10,22 


Q2359-02 


2.1539 


20.30±0.10 


1 


-1.53±0.10 


1 


-1.83±0.10 


1 


-1.53±0.10 


10,22 


Q0000-2619 


3.3901 


21.41±0.08 


5 


-1.68±0.20 


1 


-2.11 ±0.09 


1 


-1.68±0.20 


3,10,12,20 


PSS0007+2417 


3.7045 


20.55±0.15 


2 


> -1.69±0.15 


4 


-1.57±0.26 


3 


-1.27±0.27 


34 


PSS0007+2417 


3.4960 


21.10±0.10 


1 


-1.52±0.11 


4 


-I.73i0.il 


1 


-1.52±0.11 


34 


PSS0007+2417 


3.8382 


20.85±0.15 


6 


-2.12±0.28 


1 


-2.39±0.15 


4 


-2.12±0.24 


34 


Q00 10-002 


2.0250 


20.80±0.10 







1 


-1.28±0.11 


2 


-1.16±0.12 


31 


J0013+1358 


3.2814 


21.55±0.15 


1 


-2.06±0.16 


1 


-2.66±0.15 


1 


-2.06±0.16 


53 


Q00 13-004 


1.9730 


20.83±0.07 


4 


-0.70±0.08 


1 


-1.47±0.08 


1 


-0.70±0.08 


27,31 


BR0019-15 


3.4389 


20.92±0.10 


1 


-1.01±0.11 


4 


-1.53±0.11 


1 


-I.01zb0.ll 


10,22 



Note. — This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding 
its form and content. 

a = No measurement; 1= Si measurement; 2 = Si lower limit; 3 = Si upper limit; 4 = S measurement, 5 = measurement, 6 = combined alpha 
limit. 

k = No measurement; 1= Fe measurement; 2 = Fe lower limit; 3 = Fe upper limit, 4 = combined [Fe/H] limits; 4 = [Ni/H]— 0.1 dex; 5 = 
[Cr/H]-0.2 dex; 6 = [Al/H], 7 = combined Fe limit 

c 1 = [a/H]; 2=[ Fe/H]. with [M/H] = [F e/H]+0.3 dex; 3 = co mb ined [a/H] limit s, 4 = combined [Fe/H] l i mits 

d References: 1 Wolfe et al. (1994). 2Meyeretal. (1995). 3 Lu et al. (1996). 4 Prochaska & Wolfe 1 1996), 5 Prochaska & Wolfe :1997a), 
6Boisseetal. i 1998), 7 Lu et al. (1998), 8 Lopez et al. i 1999), 9 Pettini et al. : 1999), 10 Prochaska & Wolfe 1 1999). 1 1 Churchill et al. 
i 2000), 12 Molaroetal. (2000), 13 Petitjean et al. (2000), 14 Pettini et al. (2000), 1 5 Prochaska & Wolf e (2000). 1 6 Rao & Turnshek 
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and here only include the 207 DLAs with z > 1.5, as z < 1.5 
DLAs are selected for the presence of strong metal lines and 
are therefore biased. In Figure [8] we plot the metallicity dis- 
tribution function for DLAs, and the black line is a Gaus- 
sian fit with a mean [M/H] of -1.51 and a width a = 0.57. 
We note that this dispersion is not due to the uncertainties 
in the measurements, with a mean uncertainty of 0.13±0.05. 
In addition, we consider the distribution of a/Fe values for 
DLAs with z > 1 .5 that have an abundance of both an alpha 
and an iron peak element. For these DLAs, we also require 
that [a/H] < -1, a s Fe is depleted by dust in DLAs with 
higher metallicities ([Vladilo et al .1 120 1 lb . Figure [9] plots the 
[a/Fe] distribution function for 154 DLAs, and the black line 



is a Gaussian fit with a mean [a/Fe] of 0.30 and a width 
a = 0.16. The mean uncertainty of [a/Fe] is 0.08±0.05, 
which is smaller than the dispersion of [a/Fe], suggesting 
there may be intrinsic variations between DLAs. Nonethe- 
less, we also calculate the standard deviation of the mea n, and 
find that the mean value of [a/Fe]=0.30±0.02 (see also CHI ). 

3.5. Comparisons of Elemental Abundances 

We use a number of different ions for determining our el- 
emental ab unda nces in the order O, S, Si, Zn, and Fe as de- 
scribed in 33.21 In this subsection we justify our practice of 
mixing these elements in our analysis of DLA metallicities 
by comparing S, Si, and Fe to Zn in Figure [10] These mea- 
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Metallicity [M/H] 

Figure 8. Metallicity distribution function for 207 DLA metal abundances, 
[M/H], including abundances from both the literature and our new observa- 
tions. Only DLAs with z, > 1-5 are used since z < 1.5 DLAs are metallicity 
biased. The black line is a Gaussian fit to the data, with a mean [M/H] of - 1 .5 1 
and a width a = 0.57. Note the very low probability of detecting a system at 
[M/H] < —3 if this Gaussian is a proper description of the gas metallicity at 
low values. 
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Figure 9. Alpha enhancement, a/Fe, distribution function of 154 DLAs, in- 
cluding abundances from both the literature and our new observations. Only 
DLAs with z > 1 .5 and [M/H]< -1 are used since z < 1 .5 DLAs are metallic- 
ity biased, and Fe measurements with [M/H]> —1 are depleted by dust. The 
black line is a Gaussian fit to the data, with a mean a/Fe of 0.30 and a width 
ct = 0.16. 

surements are b ased on a different set of DLAs than those de- 
scribed in 33.41 These measurements are based only on very 
high resolu tion measurements from HIRES or UVES as com- 
piled by IVladilo et alj ( 1201 11). with additio nal data from the 
present study and from Wolfe et al.| (120081) . We note that for 
these comparisons we are not concerned about any metallicity 
biases, and therefore use all data that is available. 

In Figure [10] we compare different elemental abundances, 
and any any observed trends may be related to dust, nucle- 
osynthesis, and even ionization. Panel (a) in Figure [TOl plots 
[S/Zn] versus [Zn/H], Panel (b) plots [Si/Zn] versus [Zn/H], 
Panel (c) plots [Fe/Zn] versus [Zn/H], and Panel (d) plots 
[Si/S] versus [Si/H]. A standard Kendall Tau test shows that 
they are all correlated, with the strongest anti-correlation in 
Panel (c) for [Fe/Zn] and [Zn/H]. However, we must correct 
the correlation coefficients for the fact that the ratios com- 
prising the x and y axes contain the same elements. For 
instance, Panel (a) in Figure [10] contains Zn in both axes, 
and therefore the Kendall Tau test should find a correlation 
even if there is none. To account for this effect, we calcu- 
late partial correlation coefficients following the prescription 



in lAkritas & Siebertl (119961) . and find that only [Si/Zn] ver- 
sus [Zn/H] and [Fe/Zn] versus [Zn/H] (Panels (b) and (c)) are 
anti-correlated at greater than 3 a, However, we caution the 
reader that for all but [Fe/Zn] versus [Zn/H], we are looking 
for trends close to the limits of what the data can provide. 
While we correct for the partial correlation, a small number 
of outliers can effect the correlation. 

First, we find a minor 2-3 a anti-correlation in Panel (a) for 
[S/Zn] and [Zn/H]. While this is not very significant, there 
does appear to be a set of values above the dashed line at low 
metallicity, and below the dashed line at high metallicity. This 
may be due to a nucleosynthetic effect, or due to small num- 
ber statistics. Panel (b) shows what appears to be an anti- 
correlation between [Si/Zn] and [Zn/H], and to investigate 
this we therefore also compare Si with S in panel (d). Panel 
(d) shows a 2cr anti-correlation between [Si/S] and [S/H], al- 
though this anti-correlation is sufficiently weak and near zero 
that there may be no trend. In fact, other than a few outliers, 
there is no obvious correlation visible by eye. This suggests 
that the minor anti-correlation observed in Panel (b) may be a 
selection effect, or some currently not understood effect. We 
therefore conclude that if there is any dust depletion in Si like 
in the local ISM, it is a small effect, otherwise it would have 
been clearly observed in panel (d). We note that it is of course 
possible that any individual DLA is depleted, which may ex- 
plain the two outliers in the figure. Lastly, [Fe/Zn] versus 
[Zn/H] are anti-correlated as expected because Fe is depleted, 
and is expected to show higher depletion at higher metallic - 
ities. In fact, at low metallicities, Panel (c) shows that there 
is a flattening of the [Fe/Zn] trend, suggesting that depletion 
in Fe is not large for low metallicities, as found previously by 
IVladilo et all feoi ll) . We therefore use S, Si, Zn, and Fe to 
measure elemental abundances, although only use Fe at lower 
metallicities. 

When using different elements for determining metal- 
licities, a correction must be applied when combining a- 
enhanced elements with Fe-peak elements (see §4.4). While 
S and Si are a-enhanced in the Milky Way, the nature of 
Zn is less clear. The comparisons in Figure [10] shows that 
[Fe/Zn] is offset from even at low metallicites, while [Si/Zn] 
and [S/Zn] are around 0, suggesting that Zn is enhanced in 
a way similar to a-elements, rather than Fe-peak elements. 
This supports treating Zn as an a-element in our metallic- 
ity evolution plots, as well as in our measurement of the a- 
enhancement of DLAs. 

We can also use this same high resolution dataset as 
in Figure [10] to explore the a-enhancement of DLAs with 
higher precision than in Figure [9] Figure QT| plots the a- 
enhancement, [a/Fe] as a function of the metallicity using Si, 
S, and Zn as alpha elements. This shows that while there is a 
large scatter in the [a/Fe] values, at low metallicities ([a/H] 
< -1 .0) the [a/Fe] distribution flattens out to a mean value of 
0.27 ± 0.02, where the uncertainty is the standard deviation 
of the mean. A Gaussian fit to [a/Fe] has a mean of [a/Fe] 
of 0.26 and a width a = 0.12, which matches that found for 
the larger lower resolution unbiased sample used in most of 
this paper. At higher metallicities, there is a correlation of 
[a/Fe] with [a/H], possibly caused by the depletion of Fe, 
which is why we limit ourselves to lower metallicity DLAs 
([M/ H] < -1) when comparing the a-enhancement of DLAs 
(see {33.71 i. This systematic offset of [a/Fe] from in the 
regime where Fe is not depleted suggests that DLAs are al- 
pha enhanced. 
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Figure 10. Comparison of different elements used to determine metallicities of DLAs. Panel (a) compares S and Zn, Panel (b) compares Si and Zn, Panel (c) 
compares Fe and Zn, and Panel (d) compares Si and S. The Kendall Tau test shows that they are all technically correlated. The correlation in Panel (a) may be 
spurious, and that of in Panel (b) may be due to selection effects. Panel (c) is clearly correlated, although at low metallicities flattens out. Panel (d) suggests that 
Si is not more significantly depleted than S, aside from a few special cases. 




[a/H] 

Figure 11. a-enhancement, [a/Fe], as a function of the metallicity using Si, 
S, and Zn as alpha elements. At low metallicities ([a/H] < -1.0) the [a/Fe] 
distribution flattens out to a mean value of 0.27 ± 0.02. At higher metallici- 
ties, there is a correlation of [a/Fe] with [a/H], which may be caused by the 
depletion of Fe. This systematic offset of [a/Fe] from in the regime where 
Fe is not depleted (low metallicity) suggests that DLAs are alpha enhanced. 



3.6. Evolution of Cosmic Metallicity 

We trace the build-up of heavy elements in neutral gas 
across cosmic time by investigating the metallicity evo- 
lution of DLAs out to z ~ 5. The following analysis 
presents the results under the assumption and expectation 



that dust obscur ation has a modest effect on the majority of 
DLA sightlines dEllison et al.l2001b|:lHerbert-Fort et al 120061: 
Uorgenson et al.ll2006tlMeiring et al]|2006UKhare et al.ll2007l) . 
In Figure [12] we plot the metallicity as a function of redshift 
for both our measurements and those from the literature with 
a total of 241 abundances. Specifically, we plot the metal- 
licities from the literature, and overplot our new [M/H] mea- 
surements as stars, with ESI measurements as green stars and 
HIRES measurements as gold stars. Of these new [M/H] mea- 
surements, 30 are at z > 4. 

There are three interesting features in this plot: 1) there 
exists a large dispersion in [M/H] at all redshifts, 2) there ap- 
pears to be a metallicity "floor" at [M/H]« -3, and 3) the av- 
erage [M/H] decreases with increasing redshift. The large dis- 
persion in Figure[T2]is not an observational error, but is due to 
intrinsic scatter amongst the DLAs. The dispersion in [M/H] 
is ~ 0.5 dex, and does not appear to evolve with redshift. The 
metallicity "floor" occurs at [M/H]« -3, with no DLAs exist- 
ing at [M/H] < -2.8 despite the sensitivity of our spectra for 
finding DLAs with [M/H] lower than -4. The implicatio ns of 
the scatter and the metallicity"floor" will be discussed in jj4.ll 

The evolution of [M/H] is investigated by computing the 
cosmological mean metallicity, (Z), where 



(Z): 



log 10 [M/ " ] W(HI),-/ 5>(HI) ; ] , (1) 



14 Rafelski, Wolfe, Prochaska, Neeleman, Mendez 

OGyr 7.7Gyr 10.3Gyr 11.4Gyr 12Gyr 12.3Gyr 



Literature + 




12 3 4 

Redshift (z) 



Figure 12. DLA metal abundance, [M/H], versus redshift. Red plus signs from the literature. The metallicities from ESI are green stars and those from HIRES 
are gold stars. The 9 blue crosses with error bars show cosmic metallicity, (Z), where horizontal error bars are determined such that there are equal numbers of 
data points per redshift bin. The vertical error bars represent 1<t confidence levels given from our bootstrap analysis. Black straight line is linear fit to the (Z) 
data points. 



and i represents each bin of DLAs as a func tion of redshift 
(iLanzetta et al.l[l99l IProchaska et alJl2003bl) . Stated differ- 
ently, (Z) = log(QM/fiff/)-log(OM/Oi//)0, where fi M and 
f!,Hi are the comoving densities of metals and of atomic hy- 
drogen. We bin the data such that each bin contains an equal 
number of DLAs (26 DLAs/bin), except for the lowest and 
highest redshift bins which contain 34 and 25 DLAs respec- 
tively. We maximize the number of bins while keeping the 
resultant values uncorrelated and require a sufficient number 
of DLAs per bin such that we obtain Gaussian distributions in 
the bootstrap analysis for determining our uncertainties. (Z) 
is plotted as blue crosses in Figure[T2] and the error bars rep- 
resent 1 a confidence limits based on a bootstrap error analysis 
described below. 

The (Z) statistic is dominated by the systems with the 
largest Nm and [M/H] values (although we find no correla- 
tion of A^hi with [M/H]), and therefore the measurement un- 
certainty is determined by sample variance rather than sta- 
tistical error. We calculate the uncertainties on (Z) by us- 
ing a bootstrap method, where we replace a random fraction 
(1/e w 37%) of the DLAs in each bin with other DLAs from 
that bin. We repeat this 1000 times, and then take the stan- 
dard deviation of all the column density weighted metallic- 
ity measurements. We chose the number 1000 because it is 
sufficiently large such that the distribution of (Z) resembles 
a Gaussian distribution, and increasing this number does not 
change the result. This method is conservative, and includes 
all the uncertainties associated with the variance of the sam- 
ple, if the distribution of the bootstrap values are Gaussian. 
We note that we have sufficiently high Nm DLAs in our sam- 
ple such that the distribution is close to Gaussian, and the re- 
sult s are insensitive to the addition of any single DLA (see 
also IProchaska et aT]|2003bl) . The only exception to this is the 
inclusion of J 1201+21 17 at z a b s = 3.79, which has an extreme 



velocity width showing absorption lines over a 1500 km s 
range. Regardless of how we bin our data, the cosmic metal- 
licity bin including this DLA is slightly elevated as evidenced 
in the z ~ 4 bin in Figure Q~2] The effect of this DLA on the 
binned cosmic metallicity is included in the bootstrap analy- 
sis, resulting in a slightly larger uncertainty for this value. 

The black dotted line in Figure [12] is a linear fit to the (Z) 
values and their bootstrap uncertainties. In order to properly 
treat the uncertainties, we perform a chi-square minimization 
fit to a logarithmic function on the linear values of the cos- 
mic metallicities and their associated uncertainties from the 
bootstrap. We find that (Z) = (-0.22±0.03)z- (0.65±0.09), 
excluding rounding errors. This is a reasonable, although 
non-physical, fit to the data to describe the slope of the evo- 
lution. This yields an 6.7 a detection in the evolution of (Z), 
which is a si gnificant improvement over the previous ~ 3a 
detections bvlProchaska et alj d2003b1).lKulkarni et alj (120051) . 
iKulkarni et al. (20071). and lKulkarni et alJ j2010h. While the 
values of t he fit are extremely s imilar to what was found pre- 
viously by Kul karni et alj (120101) . the significance of the evo- 
lution is improved upon by the larger lever arm obtained by 
sampling DLAs at higher redshift, more DLAs at lower red- 
shift, and a larger sample size. We note that the sample com- 
piled from the literat ure (see T able |3} is somewhat different 
than that compiled by Kulkarn i et alj (120101) . since we exclude 
some studies due to potential metallicity bias, which may also 
affect the results. Regardless, there is good agreement of the 
slope and intercept of these semi-independent studies. 

The larger the number of cosmic metallicity measurements, 
the better we can determine the evolution of the cosmic metal- 
licity. However, if we reduce the number of DLAs per bin too 
significantly, then our binned values will be plagued by small 
number statistics, and the bootstrap may not result in Gaus- 
sian distributions. As a test to investigate our sensitivity to 
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binning, we apply a running bin with 30 DLAs per bin. Each 
bin adds and removes one DLA while increasing in redshift, 
resulting in 202 bins. This methodology results in correlated 
values, as each DLA is used many times, and thus is not used 
as our primary technique. However, it is a good test of the ef- 
fect of binning on our data. This yields (Z) = (-0.19)z-(0.72), 
where we omit the error bars due to the heavy correlation of 
the data. The slope and intercept are very similar to the above 
values, and are well within the uncertainties of the uncorre- 
cted data. This suggests that binning is not strongly affecting 
our results. 

In the fit of the cosmic metallicity, we make no correc- 
tions for the pote ntial metallicity bias in our lowest redshift 
bin mentioned in 33.41 Specifically, the DLAs at z < 1.5 are 
generally selected based on their strong Mg II, rather than on 
the damped Ly-a line. The two selection methods are not 
equivalent, and DLAs selected for strong Mg II absorption 
syste ms are generally more metal rich than Ly-q se lected sys- 
tems dNestor et al.ll2003ll2008tlMurphv et al 120071) . Correct- 
ing for this bias of ~ 0.1 dex dNestor et al.l 120081) results in 
(Z) = (-0.19±0.03)z-(0.74±0.09), which is a 6.3a detec- 
tion in the evolution of (Z). On the other hand, removing the 
lowest redshift bin completely, and thereby only using DLAs 
with z > 1.5, results in (Z) = (-0.22±0.05)z- (0.66±0.15), 
which is a 4.6a detection in the evolution of (Z). We note 
that all the fit values ar e completely consistent w ith the pre- 
vious measurements bv lProchaska et all (1200 3b). confirming 
that the mean metallicity of the universe in neutral gas is dou- 
bling about every billion years at z ~ 3. 

3.7. Comparison of Metallicities of DLAs and Halo Stars 

DLAs are widely believed to be the progenitors of to- 
day's disk galaxies and act as neutral gas reservoirs for 
star formation at high redshifts (Nagamin e et al.l 120 04b a; 
IWolfe & Chenl l2006: Wol fe et alJl2008t iRafelski et alJl201 ll) . 
Stars forming in situ out of low metallicity DLA gas would 
result in low metallicity stars, which should be observable in 
the Milky Way today. However, previous studies found that 
the metallicity distribution of DLAs disagree with known stel- 
lar po pulations in the Galaxy dPettini et al.l 1 997t iPettini 2004, 
2006). On the other hand, recent studies of metal-poor DLAs 
suggests that their [C/O] distribution is in good agreement 
with m etal-poor halo stars (Penpr ase et al.ll2010l ICooke et al.l 
12011b . 

Here we investigate whether the increase in sample size and 
redshift yields a metallicity distribution in better agreement 
with Galactic stellar populations. We compare the metallicity 
distributions of DLAs with the thin disk stars, thick disk stars, 
and halo stars, although we mainly focus on the halo stars as 
the metallicities of the other two are systematically too high 
(as we show below). We limit ourselves to DLAs at z > 2, 
corresponding to an age of >10 billion years, similar to the 
ages of halo stars in the Milky Way. However, the results are 
insensitive to including DLAs at z > 1.5. 

Halo stars generally have low metallicities and 
are often selected ba s ed o n their metallicities (e.g. 
iBeers & Sommer-Larsenl 119951) . However, there is an 
overlap in the metallicity of different stellar component s 
in the gal axy (|Unavane et al.l 119961: IChiba & Beers! 2000). 
Therefore, IVenn et al.l (120041) select stars based purely on 
their kinematics, compile a large dataset from a number of 
publications, and provide probabilities of each star consist- 
ing of a thin disk, thick disk, or halo st ar, based on their 
kinematics. Here, we select stars from the IVenn et al.l (12004b 



sample of local stars with an 80% or better probability of 
being identified as each type of star and having both a and 
Fe element abundances measured. This results in a sample 
of 20 1 thin disk stars , 109 thick disk stars, and 207 halo 
stars (IVenn et alj|2004j: lEdyardsson et alJll993t lHanson et al.1 
1998: lFulbrightll2000l.l2002HReddy et alJ2003l) . 

While [Venn et al.l (12004) sample does not select stars based 
on their metallicity, the sample is still metallicity biased as the 
studies it uses preferentially target metal poor objects. This is 
a well known bias, and its effects a re discussed in depth for 
a different sample of halo stars by (ISchorck et al.l 120091) . In 
that case, the metallicity distribution function was offset by 
0.3 dex lower than the intrinsic distribution. While larger less 
biased samples of stellar metallicities and clas sifications ex- 
ist, th ese data have not been made public (e.g. ICarollo et al.l 
2010). In addition, while our DLA sample is not metallic- 
ity biased per redshift bin, it is slightly biased as a whole 
due to the metallicity evolution with redshift. If we had tar- 
geted additional high or low redshift DLAs, it would slightly 
change our metallicity distribution function by decreasing or 
increasing the mean metallicity respectively. Since we are rel- 
atively evenly sampled at z > 2, we expect this to be a small 
effect. Moreover, a comparison of stellar metallicities with 
those of DLAs should be approached with caution, as all the 
stars from any given stellar population in the Milky Way could 
have formed out of just a small number of DLAs. Also, while 
it is logical to assume that the Milky Way is a representative 
object, there will be variations between galaxies that affect 
the comparison. Hence, while we expect the metallicities of 
stars formed out of DLA gas to be similar to those of stars in 
the Milky Way, they do not need to be exact matches. It is, 
however, illustrative to compare the metallicity of DLAs with 
those of different stellar populations. 

In Figure [13] we plot histograms comparing the abun- 
dances of 195 DLAs at redshifts z > 2 (red), with a) thin 
disk stars, b) thick disk stars, and c) halo stars (blue). For 
the DLAs, we use the same metal abundances as used in Ta- 
ble [2] and Figure Q~2] F or the stars , we u se the a-element 
abundances tabulated in IVenn et all (120041) . The histograms 
are normalized to have equal area and a maximum value 
of 1.0. We compare the stars to DLAs with z > 2, as the 
majori ty of the halo stars are believed to have formed by 
then (Freeman & Bland-Hawthorn 2002; Bullock & Johnston 
120051; [Robertso n et al.ll2005t iJohnston et al. 2008). This cut 
includes the vast majority of our sample, and therefore has 
a negligible effect on the comparisons. A visual compari- 
son of histograms immediately rules out the possibility of any 
agreement of the metallicity distributions of the thin disk or 
thick disk stars with DLAs. However, the halo star metallicity 
distributions look similar, and the median [M/H] of the two 
distributions are basically the same, with values of -1.54 and 
-1.53 for the DLAs and halo stars respectively. 

We apply the Kolmogorov-Smirnov (K-S) test to find the 
probability of each of the distributions being drawn from the 
same parent population as the DLA metallicities. Similar to 
the visual comparison, the thin disk and thick disk stars are 
ruled out at ~ 0% probability, i.e., the null hypothesis can be 
rejected at a high confidence level. However, the halo stars 
have a probability of 3%, i.e., the null hypothesis cannot be 
rejected with more than 2a confidence. The two samples are 
therefore consistent with being drawn from the same parent 
population. Furthermore, we reiterate the warning mentioned 
above, that the two distributions do not need to match exactly. 
For instance, the metal poor tail shown here for halo stars is 
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not evident in th e metallicity d istribution function from SDSS 
(see Figure 2 in ICarollo et al.ll2Q10t) . and is due to the metal 
poor bias mentioned previously. In any case, this compari- 
son hints at an overlap between the metallicity distributions 
of high redshift DLAs and halo st ars i n the Milky Way. We 
discuss the implications of this in 34.21 

Another possible comparison of DLAs and halo stars can 
be made by studying the chemistry of the two popula- 
tions. Specifically we can compare the ratio of a-elements 
to Fe abundances to determine whether gas is a-enhanced. 
a-elements are produced in high mass stars and then are 
ejected by Type II supernovae (SNe), while Fe is produced 
in both Type II and Type la SNe. Therefore, if stars are 
formed shortly after the Type II SNe explode, but before 
the Type la SNe do, they will have enhanced [a/Fe] ratios 
(iMatteucci & Recchill2001l) . While the exact lifetimes of SNe 
is an active area of research, Type II SNe (and other a-element 
producers such as Type lb and Ic SNe) generally need 10 6 -10 7 
years t o explode, while Type la SN e need more than 10 8 -10 9 
years (Kobayashi & Nomoto 2009). What is important is that 
Type la SNe take significantly more time to produce the Fe 
elements than the other SNe require to produce a-elements. 

In Figure[l4] we plot histograms of the [a/Fe] ratios for the 
207 halo stars (blue) and a) 138 z > 2 DLAs that have abun- 
dances for both an Fe and an a-element, and b) a subset 115 
of those 138 DLAs that also have [M/H] < -1. We note that 
it is possible that the Fe abundances are depleted onto dust 
grains, especially for the higher metallicity systems. How- 
ever, at lower metallicities, the [q /Fe] ratio of DLAs is rela- 
tively constant (s ee Figure QTJ and Prochaska & Wolfe! 120021 
Wolfe et al.l2005h . If there was dust depletion in the sample, it 
would move the red histogram to the right in Figure [14k The 
sample in Figure [14b reduces such depletion issues, as the 
lower metallicity DLAs will have Fe less or not at all depleted 
(see FigurefTTll. In comparison, the mean alpha enhancements 
of the thin disk and thick disk are ~ 0.16 and ~ 0.20 respec- 
tively, inconsistent with those of DLAs. 

We apply the K-S test to the [a/Fe] ratios for DLAs and 
halo stars, and find the probability that the two distributions 
are drawn from the same parent population to be 3% for the 
138 DLAs in FigureUli, and 40% for the 1 15 DLAs in Fig- 
ure fl4b . We therefore cannot reject the null hypothesis with 
more than la confidence, and find that the [a/Fe] ratio dis- 
tribution for the DLA sample and the halo stars are com- 
patible. This suggests that both the halo stars and the DLA 
gas are a-enhanced by rj 0.3 dex, suggesting that halo stars 
could be formed out of D LA gas. This result is different than 
some previous studies (e.g.lPettini et al.|1999HCenturi6n et all 
20051 iMolaro et al Jl200(ilLedoux et al.ll2002at I Vladiloll2002t 
Nissen et al JI2004L 120071) mainly due to a differ ent interpreta- 
tion of the behavior of Zn (see $3.21 and §13. 5K and previous 
small number statistics. For instance. iNissen et al.l (120041) find 
no [S/Zn] enhancement in DLAs, and therefore conclude that 
there is no a-enhancement. However, we find that Zn behaves 
like an a-element, and therefore does not constrain the [a/Fe] 
ratio. 

4. DISCUSSION 

4.1. Metallicity Evolution 

The principal result of this study is that we find a continued 
decrease in the metallicity of DLAs with increasing redshift 
out to z ~ 5. We improve the significance of this trend to 6.7a 
and extend it to higher redshifts. These observations of the 
metallicity evolution tightly constrain the star formation his- 
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Figure 13. Histograms comparing abundances of 161 DLAs at z > 2 with a) 
201 t hin disk stars, b) 109 thick disk stars, and c) 207 halo stars JVenn et all 
2004). The histo grams are normalized to have equal area and a maximum 
value of 1 .0 The red left shaded region represent the DLAs and the blue right 
shaded region represent the thin disk, thick disk, and halo stars. The only 
population of stars in the Milk Way that is consistent with being drawn from 
the same parent population are the halo stars. 

tory of z ~ 5 galaxies and the processes that transport metals 
from star-forming regions to the ambient ISM. 
Figure [l5]compares current models to our measured abun- 
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(b) z > 2 and [M/H] < -1 DLAs 

Figure 14. Histogram comparing a/Fe ratios of 207 halo stars with a) 138 
DLAs at z > 2 and b) 115 DLAs at z > 2 with [M/H] < -1. The DLAs are 
selected such that they have both a and Fe elemental abundances measured. 
The histograms are normalized to have equal area and a maximum value of 
1 .0. The red left shaded region represents the DLAs and the blue right shaded 
region represents the halo stars. The cn/Fe ratios of DLAs is consistent with 
that of halo stars. 



dance evolution. We show the same data and best fit cos- 
mic metallicity as in Figure Q~2] and over-plot three mod- 
els from simulations. The red diamonds are mean metallici- 
ties of DLAs extracted from th e 48 Mpc/h momen tum-driven 
wind simulation present ed in iDave et al.l (1201 11) . using the 
self-shielding criterion in Popping et al. (2009), in cluding an 
molecular fraction calculation based on Blitz & Rosolowskv 
The 



(2006). The green triangles are from Fumagal li et al.l 
( 201 lbl) based on simulations by iCeverino et al.l (1201 Oh us- 
ing an atomic fractio n calculation based on the formalism by 
Krumholz et al. (2009). Lastly, the gold squares are directly 
from|Cen|(|20T2). 

The normalization of the metal abundance of DLA gas in 
the simulations are remarkably similar to the observed val- 
ues for all three models. In addition, the evolution of the 
me tal abundance is v ery si mila r for tw o of them. The models 
by IDave et al] (1201 lh and ICenl (120121) agree within their un- 
certainties, although they strictly do have ste eper slopes than 
the data . On the other hand, the models by Fumagalli et al. 
d201 lbl) show no real evolution in the DLA metallicity, which 
may be due to sampling a small number of halos, or other 



issues in the simulation. In all, it is reassuring that differ- 
ent models with different hydro-formulations seem to broadly 
agree on the metal content, and match the measured values. 
However, although these models treat metal diffusion, mix- 
ing, and feedback differently, the broad agreement between 
the models suggests that perhaps the metals are mostly locked 
in the high density gas where stars are produced (private com- 
munication, Michele Fumagalli 2012). 

In addition to the continued decrease in metallicity, Figure 
[T2]reveals two interesting features. The first is the large dis- 
persion of [M/Hl at z > 4 that matches the results at lower 
redshift (Prochaskaet al. 2003b). This result is not an obser- 
vational uncertainty but is due to intrinsic scatter amongst the 
DLAs. The scatter does not appear to evolve with redshift, 
suggesting that it is intrinsic to the DLAs, possibly due to the 
wide range of masses of the galaxies hosting DLAs (Neele- 
man et al. 2012 in prep), as suggested b y the mass-metallicity 
correlation found bo th at low redshift flT remonti et al.l 120041) 
and at high redshift ( Erb et al.ll2006l iProchaska et al.ll2008l) . 

The second feature is a metallicity "floor" at [M/H]«-3, 
which is also consistent with the results at lower redshift 
(IProchaska et alJl2003bt iPenprase et all 120101) . We note that 
we cannot currently rule out the alternative possibility that 
rather than a floor, this cutoff is due to the tail of the presum- 
ably gaussian distribution not bein g adequately sample d (see 
Figure [HJ. This is investigated in Penpr ase et al.l (120101) . who 
find that if it is the gaussian tail, the probability of finding a 
DLA with [M/H]< -3 is extremely small. 

If the observed metallicity "floor" is real, then it is likely a 
physical lower limit of the metallicity of DLAs out to at least 
z«5, which may continue out to even larger redshifts. We 
note that if there is gas with very low metallicity, it has a very 
small cross-section at Nm> 2 x 10 20 cm" 2 . The observed DLA 
metallicities are systematically higher t han those of the Ly - 
a forest clouds which have [M/H]=-3 ( Ag uirre et al.1 12004). 
and therefore even the lowest metallicity DLA gas is part of 
a distinctly different population. Moreover, the extension of 
this lower limit of [M/H] at higher redshift has implications 
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Figure 15. Comparison of the DLA metal abundance evolution from simu- 
lations to the data. The gray plus signs are all DLAs used in this study, the 
blue crosses with error bars show the cosmic metallicity, (Z) , and the black 
straight line is a linear fit to the (Z) data points. The red diamonds are mean 
metallicities of DLAs extracted from the sim ulation presented inlDave et all 

12011) . The green triangles by models from Fumagalli et al. (2011b) based 
on sim ulations by |Ceverino et al. 1 20H|). The gold squares are from ICenl 

12012) . All the error bars are la confidence levels. The normalization of 
the metal abundances all agree with the data, while the slope in the evolution 
of the metal abundance is somewhat different, although consistent within the 
uncertainties. 
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for the primordial gas in these galaxies. Either the primor- 
dial gas does not exist in the neutral phase in high redshift 
galaxies, or always coin cides with a metal enriched region 
dProchaska et al.l l2003bh . Simulations of primordial metal 
enrichment including Population III stars show that a sin- 
gle pair-instability supernova can enrich a halo to [M/H]ss-3 
dWise et ail 120121) . which would explain and be consistent 
with the observed "floor". Such a pre-enrichment of the ISM 
by Popullation III stars resulting in a minimum metallicity 
would also improve the a greement of simul ated galaxies and 
low-mass dwarf galaxies (Tassis et al. 2 0121) . 

While the evolution of (Z) is linear in redshift, it is non- 
linear in time, with the slope, d(Z)/dt, becoming significantly 
steeper at earlier cosmic times. In Figure [16] we plot the 
metallicity as a function of time rather than redshift. We trun- 
cate the plot at z = 1 .6 in order to focus on the higher red- 
shift part of the plot which contains the most rapidly chang- 
ing metal abundances. The black dotted line is the linear fit to 
(Z) versus redshift, shown as a function of time, and clearly 
shows the steepening of (Z) going back in time. While this 
trend line is not a physical model, it shows the importance of 
obtaining metallicities of the highest redshift DLA systems, as 
undertaken in this study. Moreover, we show the metallicity 
prediction from sim ulations of prim ordial metal enrichment 
as a black diamond (Wise et al. 2012), which would require a 
drop in metallicity below what is extrapolated from the data. 

An intriguing property of the [M/H] values in Figure [T2"lis 
that the highest redshift DLAs in our sample appear to have 
metallicities below the linear trend in (Z). Specifically, 6 of 
the 7 DLAs at z > 4.7 have [M/H] values below the linear fit 
((Z)= -1.6 at z = 4.8), with a column density weighted mean 

(Z) = -2.1+0 2 ( The DLA that does not is the l° west redshift of 
the group). We note that the linear fit is not a physical model, 
and therefore does not constrain our expectations. However, 
the physical models shown in Figure[l5]do not currently show 
a sudden drop in metallicity, so if real, this would mark a de- 
parture of these models from the data. Although the sample 
size is currently small, and this result may just be due to small 
number statistics, future work should test whether this depar- 
ture of (Z) below the linear fit persists at at z > 4.7. 

4.2. Comparison with Halo Stars 

The metallicity measurements at high redshift also enable 
us to investigate whether any known stellar populations in our 
Galaxy have the same metallicity distribution as the high red- 
shift DLAs. Previously, no known population w as known to 
have the same metallicity distribution as DLAs dPettini et al.l 
1 1 9971: lPettinill2004L 120061) . This was because the metallicity of 
thin and thick disks were too metal rich for DLAs, while the 
halo stars were too metal poor. However, we showed in § 13.71 
that it is possible that the metallicity distribution of DLAs at 
z > 2 is compatible with that of halo stars. Specifically, we 
cannot rule out, with a high degree of confidence, that DLAs 
are the progenitors of halo stars. We stress that this is the first 
time that the metallicity distribution of any known population 
of stars have been shown to be consistent with being drawn 
from the same parent population as DLAs. 

There are two primary differences between our study and 
previous studies, which compared the D LA metallici t y dis- 
tributions to know n stellar populations dPettini et al.l 11997b 
lPettinill200l 12006). First, we use a-elements for abundance 
measurements (primarily S and Si), while the other studies 
use Zn. Second, we probe a higher redshift range out to z ~ 5, 
as opposed to z ~ 3.5. Both Zn and Si yield reliable abun- 



dance me asurements, which are consistent with each other 
(see §!3.51 l. although, the oscillator strengths of Zn are smaller 
than Si, making it much more difficult to detect low metallic- 
ity DLAs. Additionally, the wavelengths of Zn transitions are 
larger than Si, making Zn measurements at high redshift more 
difficult to obtain as they are redshifted into the near-infrared. 
The past comparisons of the metallicity distributions included 
upper limits of Zn together with detections, skewing the dis- 
tribution towards metallicities higher than that of the intrinsic 
population. Furthermore, by sampling higher redshifts, we in- 
clude lower metallicity DLAs, since the metallicity of DLAs 
evolve with time. Together, our sensitivity to lower metallicity 
systems and our inclusion of higher redshift systems account 
for the different metallicity distribution presented here. 

In addition to the metallicity distribution, the chemistry of 
the DLAs and halo stars is consistent, based on the compari- 
son of the [a/Fe] ratios. Both populations have enhanced a- 
elements ([a/Fe] « 0.3), suggesting that stars formed shortly 
after the ISM was enriched by type II super novae, but before 
the type la supern ovae had a chance to go off (Mat teuccil2003t 
IVenn et aT]|2004l) . Together, these results suggest that the halo 
stars formed out of gas that has undergone the same number 
of generations of star formation as DLAs, which occurred be- 
fore type la supernovae contributed significantly. It is there- 
fore possible that the halo stars in the Milky Way formed out 
of z > 2 DLA gas. 

There are two primary formation scenarios for the Galac- 
tic ha lo stars (for a more thorough review, see iHelmil 
(2008)). The first is an in-situ scenario with the formation 
of halo stars from in-falling gas towards t he Galactic center , 
which later collapse onto the plane (e.g. Egge n et ail 1 19621; 
Samland & Gerhardl 120031) . The second is an accretion sce- 
nario where a num ber of smaller galaxies merge to form the 
Galactic halo (e.g. [ Searle & Zinnlll978t iBullock & Johnston! 
120051: iMoore et al.l 12006b . The later is often considered the 
more favored scenario, due to its resemblance to the model 
of hierarchical assembly with the collapse of sm all haloes of 
cold dark matter (CDM) (IWhite & Reeslll978l) . It is likely 
that the answer lies s omewhere in-be t ween, with both pro- 
cesses playing a role (|M aiewski 1993; Chi ba & Beersl 120001: 
IHelmil 12008 1 : iFont et all 1201 11) . In addition, while Galactic 
halo stars used to be considered a single stellar population, 
recent evid ence suggests the existence of an inner and an 
outer halo (ICarollo et al.ll2007l: lLai et al J 120091: ICarollo et al.l 
2010). The formation mechanism of these two populations 
are likely different; with the inner halo formed from a com- 
bination of accretion of early rapid mergers and in-situ for- 
mation, and the outer halo formed purely from the disruption 



of lower mass satellite ga l axies and accr etion (Carollo et al 



20071: lZolotoveTaII [2009: Ca rollo et al.ll20Tot 



2010). 



Cooper et al 



The Galactic stellar halo formation scenarios are com- 
patible with halo stars formed out of DLA gas. Specif- 
ically, these scenarios include the mergers of lower mass 
systems that form the majority of the stars. While nearby 
dwarf spheroidal galaxi es were originally thought to lack 
stars with [Fe/H] ~ -3 dHelmi et alj |2006). fainter systems 
have sufficiently low met allicity that they could be contribu- 
tors to the Galactic halo dKirbv et alJl2008t iKoch et alj|2008h 
iNorris etal.ll200l iFrebel et al.ll2010t lLai et al.l 1201 11) . Such 
low metallicity stars are even now seen in more massive 
dwarf spheroidal galaxies (iKirbv et al.ll2.Ql lbl) . In addition, 
while initially the a-enha ncement of the dwarf galaxies did 
not match the stellar halo (She trone et al.ll2003l: iGeisler et al.l 
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Figure 16. DLA metal abundance [M/H] versus time similar to figure [12] The black diamond is the prediction from a hydrodynamic simulation on primordial 
metal enrichment (Wise et al. 2012). The purple square is a column density weighted mean of the 7 highest redshift DLAs at z > 4.7 in our sample, which all 
have values below the fit. Red plus signs from the literature. The metallicities from ESI are green stars and those from HIRES are gold stars. The 9 blue data 
points with error bars show cosmic metallicity, (Z), where horizontal error bars are determined such that there are equal numbers of data points per redshift bin. 
Black straight line is linear fit to (Z) data points as a function of redshift. We note that we show only the values at z. > 1 .6 to focus on the area of interest at higher 
redshift. 



2007), th ere is now agreement for low metallicity stars in such 
galaxies dFrebel et al.ll2010l [Kirby et all 1201 lab . Similar to 
dwarf galaxies that may have been halo progenitors, the a- 
enhancement and metallicity distribution of DLAs at z > 2 
match the halo stars (see Figures [T3l and [141. In fact, DLAs 
encompass both l ower mass galaxies and the outer regions 
of spiral galaxie s (Matte ucci et al.lll997t iPro chaska & Wolfe ! 
1 19981: iDessauges-Zavadskv et al.1 12004 120071: iMeiring et al] 
1201 latlFumagalli et alj 1201 lbl) and therefore a subset of 
DLAs are in fact dwarf galaxies (althou gh not all DLAs can 
be dwarf galaxies (Prochaska et al. 2005:!)). It is therefore not 
overly surprising that DLAs also share characteristics with 
halo stars, and that halo stars may have formed out of DLA 
gasQ 

4.3. Star Formation Rate Efficiency 

The metallicity evolution also has implications for the star 
formation rate (SFR) efficiency at high redshift. In general, 
the SFR per unit area is related to the gas density via the lo- 
cally established Kennicutt-Schmidt (KS) relation (Esfr oc 
E^ aj ). At z ~ 3, the in situ SFR of DLA gas was found 
to be less than 5% of w hat is expected from the KS rela- 
tion in isolated regions dWolfe & Chen! 120061). and ~ 10% 
of the expec ted rate in the outskirts of LBGs (Rafelski et al. 
120091 1201 11) . This means that the SFR of DLA gas is less 
efficient than what would be predicted by the KS relation. 
iGnedin & Kravtsovl (j2010) conclude that the main reason for 
the decreased efficiency of star formation is that the diffuse 

12 We note that if the halo stars formed out of DLA gas, then the higher 
kinematics of halo stars than DLAs would necessarily be caused by mergers, 
which randomize and increase the kinetic energy of the orbits of the stars 
formed in the quiescent DLAs. 



ISM in high redshift galaxies has a lower metallicity and dust- 
to-gas ratio, both of which are needed to cool the gas and form 
stars. In this case, since the metallicities of DLAs decrease 
with redshift, we expect that the efficiency of star formation 
should be correlated with redshift. In order to further our un- 
derstanding of the effect of metallicity on the SFR efficiency, 
and measure the evolution of the KS relation with redshift, 
measurements over a range in redshift are required. This study 
provides the metallicities of the DLA gas over a large redshift 
range, setting the stage for such studies. 

5. SUMMARY 

The purpose of this study is to measure the metallicities 
of DLAs with z > 4 in order to determine the metallicity 
evolution of DLAs out to z ~ 5. We obtained high resolution 
spectroscopy of 68 quasars using the ESI and HIRES instru- 
ments on the Keck Telescopes and find the following: 

(1) We fit A^hi Voigt profiles to 68 candidate DLAs from 
SDSS, and confirm 51 of them to have N m > 2 x 10 20 cnr 2 
(see Figures |2] and O. For the 47 z > 4 DLA candidates, 15 
are falsepositives leaving 32 confirmed DLAs with z > 4 (see 
Figure QJ, tripling the sample of z > 4 DLAs that have high 
resolution metallicity measurements to a total of 39. While at 
lower redshifts the majority of candidate DLAs from SDSS 
are confirmed, at z > 4 the number of false positives is large, 
with ~30% misidentified due to line-blending with random 
Ly-a forest clouds. Studies using DLA identifications from 
SDSS at z > 4 therefore need to apply a correction for the 
increased false positive rate for the lowest Nm DLAs before 
determining the column density distribution function, the 
total covering fraction, and the integrated mass density of H I 
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gas. 

(2) We consider the different elements available to measure 
the metallicity of DLAs, and find that the metallicity is 
determined with the highest confidence from elements in the 
following order: O, S, Si, Zn, Fe. We compare the elemental 
abundances determined from different ions in Figure [10] and 
conclude that S, Si, and Zn all result in similar metallicities. 
In addition, we find that for low metallicity DLAs, the 
metallicity can be determined from Fe, with a correction due 
to the alpha enhancement of DLAs. In Figures [9] and [TT] we 
show that DLAs are alpha enhanced, with a mean value of 
0.27 ±0.02. 



that transport metals from star-forming regions to the ambi- 
ent ISM. The 6 highest redshift DLAs in our sample all have 
both [M/H] and (Z) values below the linear fit. Future ob- 
servations should therefore focus on obtaining metallicities at 
Z > 4.7 to further constrain the metallicity evolution of DLAs, 
as they will have the largest impact. Lastly, we find that the 
metallicity distribution and a-enhancement of DLAs is sim- 
ilar to Galactic halo stars, enabling the possibility that halo 
stars form out of DLA gas. 



(3) We measure metal abundances, [M/H], for a total of 47 
DLAs (see Table [2] and Figure |5), with 30 at z > 4 (one 
of which is a proximate DLA). We find that the cosmic 
metallicity, (Z), continues to decrease with increasing 
redshift to z ~ 5 (see Figure [T2l . Specifically, we find 
that (Z) = (-0.22±0.03)z-(0.65±0.09), which is a 6.1a 
detection in the evolution of (Z) and is a sign ificant improve- 
ment ov er the previous ~ 3 a det ect ions by | Prochaska et al.l 
(2003b), Kulk arnietaLl J20051). iKulkarni et aTr i2007), 
and iKulkarni et all d2010h . If we remove our lowest 
redshift bin which may be slightly metallicity biased, 
and thereby only using DLAs with z > 1.5, we find that 
(Z) = (-0.22±0.05)z-(0.66±0.15), which is a4.6er detection 
in the evolution of (Z). The fit values are completely con- 
sistent with the previous measurements by Prochaska et al. 
(2003b), confirming that the mean metallicity of the universe 
in neutral gas is doubling about every billion years at z ~ 3. 
Additionally, the normalization and evolution of the metal 
abundance of simulations agree remarkably well with our 
data, as shown in Figure [TBI 

(4) We find that the large dispersion of [M/H ] of ~ 0.5 dex 
measured at z < 4 (Prochas ka et al.l 12003b!) continues out 
to z ~ 5 (see Figure [6j. This result is not an observational 
uncertainty, but is due to intrinsic scatter amongst the DLAs. 
The scatter does not appear to evolve with redshift, and could 
be due to the wide range of masses of the galaxies hosting 
DLAs (Neeleman et al. 2012 in prep). 

(5) We find that t he metallicity "floor" at [M/H]«-3 
measured at z < 4 (Prochas ka et al.l 12003b!) continues out 
to z ~ 5 (see Figure |6). Stated differently, we find no 
DLAs with [M/H] < -2.8 despite the sensitivity of our 
spectra for finding DLAs with [M/H] > -4.5. This is 
likely a physical lower limit of the metallicity of DLAs out to 
at least z ~ 5, which may continue out to even larger redshifts. 

(6) We find that the metallicity distribution and the a/Fe 
ratios of z > 2 DLAs is consistent with being drawn from 
the same parent population as halo stars (see Figures [13] and 
IT4t . This is the first time that the metallicity distribution of 
any known stellar population has been shown to be consistent 
with being drawn from the same parent population as DLAs. 
It is therefore possible that the halo stars in the Milky Way 
formed out of gas that commonly exhibits DLA absorption at 
Z>2. 

Altogether, we have measured the evolution of the metal 
abundances of DLAs out to z ~ 5, which tightly constrain 
the star formation history of z ~ 5 galaxies and the processes 
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A. ONLINE MATERIAL 

Table |4] provides the column density measurements for the 
entire new DLA sample. Column (1) gives the quasar coordi- 
nate name obtained from the SDSS survey, column (2) gives 
the DLA absorption redshifts determined from the metal tran- 
sition lines. Column (3) gives the Ion measured, column (4) 
gives the wavelength of said ion, column (5) gives the oscil- 
lator strength used, and column (6) gives the instrument used. 
Column (7) gives the velocity interval over which the equiva- 
lent width and column density are measured, and column (8) 
gives the rest frame equivalent width. Column (9) gives the 
measured column density, and column (10) gives the adopted 
column density value combining measurements fo multiple 
transitions. 

Figure set 17 shows the velocity profiles of element transi- 
tions for the DLAs in Table|4] Every column density measure- 
ment in this table is shown in these Figures. These velocity 
profiles are a mix of HIRES and ESI data, as specified in Table 
|U although it should be quite obvious based on the resolution 
of the data in the Figures. 

Fig. Set 17. Velocity profiles of DLAs. 
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Table 4 

Ionic Column Densities 



QSO 
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Ion 
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log/ 
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(km s _1 ) 
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(mA) 


log N 


log Nadopt 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


J0040-0915 


4.7394 






















CII 


1334.5323 


-0.8935 


ESI 


[-80, 100] 


449.4 ±8.0 


> 14.64 


> 14.64 






CIV 


1548.1950 


-0.7194 


HIRES 


[-40,70] 


126.9 ±21.9 


> 13.60 


13.77 ±0.03 








1550.7700 


-1.0213 


HIRES 


[-40.70] 


99.7 ±7.5 


13.77 ±0.03 








OI 


1302.1685 


-1.3110 


HIRES 


[-40. 70] 


258.5 ±7.5 


> 14.91 


> 14.91 






A1II 


1670.7874 


0.2742 


HIRES 


[-40,70] 


189.7 ±18.1 


> 12.97 


> 12.97 






Sill 


1304.3702 


-1.0269 


HIRES 


[-40,70] 


129.4 ±7.8 


> 14.16 


> 14.18 








1526.7066 


-0.8962 


HIRES 


[-40,70] 


229.5 ±22.1 


> 14.18 








Si IV 


1393.7550 


-0.2774 


HIRES 


[-40,70] 


182.7 ±9.8 


13.49 ±0.05 


13.34 ±0.04 








1402.7700 


-0.5817 


HIRES 


[-40.70] 


51.4±10.4 


13.19±0.08 








Fell 


1608.4511 


-1.2366 


HIRES 


[-30,70] 


103.0±13.1 


14.05 ±0.06 


14.05 ±0.06 








1611.2005 


-2.8665 


ESI 
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Note. — This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding 
its form and content. 
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Figure 17. Velocity profiles of metal transitions for each DLA in Table[TJ These figures are available in the online Figure set. 
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